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Abstract 

We determine the Tevatron's reach in supersymmetric parameter space in trilepton, like- 
sign dilepton, and dilepton plus tau-jet channels. We critically study the standard model 
background processes. We find larger backgrounds and, hence, significantly smaller reach 
regions than recent analyses. We identify the major cause of the background discrepancy. 
We improve signal-to-noise by introducing an invariant mass cut which takes advantage 
of a sharp edge in the signal dilepton invariant mass distribution. Also, we independently 
vary the cuts at each point in SUSY parameter space to determine the set which yields 
the maximal reach. We find that this cut optimization can significantly enhance the 
Tevatron reach. 
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1 Introduction 



For at least the next 6 years the Fermilab Tevatron will remain the highest energy collider in 
the world. The Tevatron upgrade will provide an exciting opportunity for discovering physics 
beyond the Standard Model. The hadronic environment at the Tevatron presents a number 
of challenges and extracting new physics signals can be difficult. In this respect, signatures 
with low Standard Model (especially QCD) backgrounds are extremely valuable, as they may 
provide our best opportunity for finding new physics before the LHC turns on. 

Supersymmetry (SUSY) |^ has been fascinating particle physicists for more than 25 years. 
It seems an intrinsic component of theories unifying gravity and gauge interactions such as 
string theory, M-theory or supergravity, and has played an important role in the 'second 
string revolution' of the last few years. The minimal supersymmetric extension of the Standard 
Model (MSSM) is a well-defined, renormalizable and calculable model, which offers a technical 
solution to the hierarchy problem [0], if the masses of the superpartners of the standard model 
particles are of order the weak scale. Our belief that supersymmetry might be relevant at 
energy scales accessible at present colliders is reinforced by the successful gauge coupling 
unification [P]. Also, due to decoupling the MSSM is generally in agreement with precision 
data 1^. In addition, a generic prediction of the MSSM is the existence of a light Higgs 
boson 1^, P], which is preferred by fits to data [0. In summary, the MSSM is a well-motivated 
extension of the Standard Model, which has a very rich and interesting phenomenology f^. 

Because of the relatively low (compared to the LHC or NLC) center of mass energy and 
integrated luminosity in Run II, the Tevatron is able to explore only the low end of the 
superpartner spectrum. Searches for colored superpartners (squarks and gluinos) are done in 
jetty channels, which suffer from relatively large backgrounds. On the other hand, SU(2)- 



gaugino pair production leads to a unique clean trilepton signature |T0[, which has been 
considered a 'gold-plated' mode for SUSY discovery at the Tevatron. Both CDF and DO have 
already performed Run I trilepton analyses [jll|. In light of the importance of this channel in 
Run II, it is clamant to 
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have a reliable estimate of both signal and background rates. We would prefer to deter- 
mine background rates from data, but until Run II we primarily rely on Monte Carlo 
simulations. ISAJET and PYTHIA have been two of the most commonly used event 
generators in SUSY analyses. While there is a reasonable agreement for the signal, 
PYTHIA-based studies |l^ have obtained larger values for the trilepton backgrounds 
(mainly WZ and ZZ) than ISAJET-based analyses |P, |13|, This discrepancy was 



noticed and discussed in the TeV 2000 Report [|T^, but was attributed to the different 
lepton rapidity cuts used in the various analyses. 

use an optimized set of cuts, which will maximize the Tevatron reach. A first step in 



this direction was taken in Refs. [|T4|, where softer lepton px cuts have been proposed, 



thus enhancing signal over background throughout a large part of parameter space. 

include next-to-leading order (NLO) corrections to the production cross sections. The 
corrections to diboson production [|T^, ti production [|r^ and Drell-Yan |T^], have been 
known for some time, and the corrections will soon be available for chargino-neutralino 
production as well 



identify regions of parameter space where the reach via the trilepton signature is dimin- 
ished and try to find an alternative search strategy in those regions. An example of 
this sort is the large tan/3 region with light sleptons, where one often finds that both 
the chargino and the neutralino decay predominantly to tau leptons. Then the trilepton 
signal has a very small branching ratio and the leptons are quite soft, which can make 
it unobservable at Run II, even for chargino masses as low as 100 GeV. In this case it 



is possible to recover sensitivity by considering alternative signatures with tau jets |jT2 
Another alternative to the trilepton signature is the inclusive like-sign dilepton channel 
, where the signal acceptance is increased by not requiring the odd-sign lepton in the 
event. 

In this paper, we shall try to address most of these issues. We perform detailed Monte Carlo 
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simulations of signal and background using both PYTHIA and ISAJET and explain the cause 
for the largest background discrepancy. We also determine the maximum reach by applying 
an optimal set of cuts at each point in the supersymmetric parameter space. 

Also, we make use of the presence of a sharp edge in the dilepton invariant mass distribution 
of the signal by applying a more restrictive invariant mass cut, thus reducing the WZ and 
ZZ backgrounds. As the NLO corrections to gaugino production are not yet available, we 



conservatively use leading order cross sections for all processes. Preliminary results ||T9| show 
that the /c-factor is roughly the same for both signal and background. Hence, we expect that 
the Tevatron reach will be improved once NLO corrections are incorporated. 

We show our results for the discovery potential of the upgraded Tevatron in the so-called 
minimal supergravity model (mSUGRA) p9[. This model has universal soft parameter bound- 
ary conditions at the grand unification scale, and its spectrum displays characteristic proper- 
ties. For example, the imposition of electroweak symmetry breaking results inQ |/i| > M2, so 
that the lightest chargino and lightest two neutralinos are gaugino-like. Also, the squark and 
slepton masses are generation independent, except at large tan/5 where the third generation 
masses can be lighter. This model has five input parameters: the scalar mass Mq, the gaugino 
mass Ml/2, the A-term Aq, the ratio of vacuum expectation values tan/5, and the sign of the 
fi term. We show results for /i > and Aq = 0. 

We adopt a signature driven approach by comparing and contrasting three of the cleanest 

channels for Run II - the trilepton (3L) jTO), like-sign dilepton (2L) pOf and dilepton plus 

tau jet (2L1T) [jl2| channels. The 3L channel is the long studied "gold-plated" channel. The 

2L channel has larger signal acceptance compared to 3L, but it is not a priori clear whether 

this advantage will be spoiled by the concomitant increase in the background. The 2L1T 

channel is known to be important at large tan f3, where the right-handed tau-slepton is lighter 

than the first two generation sleptons. Here we will discuss this channel at small tan (3 as well. 

We describe in detail our numerical analysis in Section |^, where we also describe the cuts 
* ^ is the Higgsino mass parameter and M2 is the soft supersymmetry breaking SU(2) gaugino mass. 
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we consider for each signature. We discuss all non-negligible backgrounds and their evaluation 
in Section ^. Then, in Section ^, we map our results for the Tevatron reach onto the parameter 
space of the mSUGRA model. We reserve Section ^ for our conclusions. 

2 Analysis 

In this section we describe our numerical analysis. We use PYTHIA v. 6.115 and ISAJET 



V. 7.42 ||22| for event generation, and the SHW v. 2.2 detector simulation package [Q, which 
mimics an average of the CDF and DO Run II detector performance. We use PYTHIA for 
the background determinations, together with TAUOLA ||2^ to account for the correct (on 
average) tau polarization in tau decays. We have made several modifications in SHW, which 
are appropriate for our purposes: 

1. We extend the tracking coverage to |?7| < 2.0, which increases the electron and muon 



acceptance, as is expected in Run II For muons with 1.5 < |?7| < 2.0, we apply the 
same fiducial efficiency as for 1.0 < \ri\ < 1.5. However, we still require that tau jets are 
reconstructed only up to \ri\ < 1.5. 

2. We retain the existing electron isolation requirement and add a muon isolation require- 
ment I < 2 GeV, where / is the total transverse energy contained in a cone of size 
AR = yA0^~+~Ar/^ = 0.4 around the muon. 

3. We increase the jet cluster Et cut to 15 GeV and correct the jet energy for muons. We 
also add a simple electron/photon rejection cut Eem/ Ehad < 10 to the jet reconstruction 
algorithm, where (Ehad) is the cluster energy from the electromagnetic (hadronic) 
calorimeter. 

4. We correct the calorimeter fix for muons. 

5. We account for an incorrect assignment of neutralino particle id's in the ISAJET trans- 
lation of STDHEP V. 4.05 [|| [|. 



^The assignment has been corrected in STDHEP v. 4.06. 



The addition of the muon isolation cut and the jet Eem/ E^ad cut allows us to uniquely resolve 
the ambiguity arising in SHW v. 2.2, when a lepton and a jet are very close. 

As we mentioned in the introduction, we show results for three of the cleanest SUSY 
channels in Run II at the Tevatron - trileptons, inclusive like-sign dileptons and dileptons plus 
a tau jet. In our analysis we consider both channel specific and channel independent cuts. In 
most of those cases, we use several alternative values for the cut on a particular variable. For 
example, we try several J^t cuts, several sets of px cuts, etc. We employ a parameter space 
dependent cut optimization: at each point in SUSY parameter space, we consider all possible 
combinations of cuts, and determine the best combination by maximizing S/ \/B. In contrast 
to superior neural network analyses, the additional CPU requirements when employing this 
simple optimization are negligible. We concede that it may not be possible to perform an 
identical analysis with real data, particularly due to trigger issues. Even so, it is useful and 
interesting to see which cuts work best in the different parts of parameter space, and to see 
how much one can gain by choosing optimal cuts. 

We first list the channel-independent cuts, which in general are designed to suppress back- 
grounds common to all three channels. 

1. Four cuts: > {15, 20, 25} GeV or no cut. 

2. Six high-end invariant mass cuts for any pair of opposite sign, same flavor leptons. The 
event is discarded if: \Mz - mi+i- \ < {10, 15} GeV; or m£+£- > {50, 60, 70, 80} GeV. 

3. Four azimuthal angle cuts on opposite sign, same flavor leptons: two cuts on the differ- 
ence of the azimuthal angle of the two highest pr leptons, \Aip\ < {2.5,2.97}, one cut 
\Aip\ < 2.5 for any ]3air leptons, and no cut. 

4. An optional jet veto (JV) on QCD jets in the event. 

We list the channel-specific pt cuts in Table |T|. In the 3L channel, the first four p^ cuts in 
the table also require a central lepton with px > 11 GeV and \ri\ < 1.0 or 1.5. 
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Channel 


Pt cuts 


3L 


Pt(^i) 






1 


11 


5 


5 


2 


11 


7 


5 


3 


11 


7 


7 


4 


11 


11 


11 


5 


20 


15 


10 



2L 






1 


1 1 


Q 


2 


11 


11 


3 


13 


13 


4 


15 


15 


5 


20 


20 



2L1T 


Pt(4) 




Pt{t) 


1 


8 


5 


10 


2 


8 


5 


15 


3 


11 


5 


10 


4 


11 


5 


15 



Table 1: Channel-specific sets of pr cuts. 



For all channels we impose a low-end invariant mass cut on any pair of opposite sign, same 
flavor leptons, m^+^- > 11 GeV. This cut is designed to suppress a number of backgrounds, 
e.g. Drell-Yan, 66, cc, and the contribution from W'^* (see below). It is common for lepton 
analyses to include a cut on the (Ar^, Ac/?) distance between any two leptons Ai? > 0.4, which 
suppresses background from 66, cc, and anomalously reconstructed cosmics. We choose not to 
include this cut, since we do not simulate those backgrounds. (Monte Carlo simulations do 
not reliably estimate the backgrounds from 66 and cc production.) We have checked, however, 
that the effect of the Ai? cut on signal and background is negligible. 

In the next section we briefly discuss the main backgrounds for the three channels. This 
will also motivate the choice of some of the cuts above. 
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3 Backgrounds 

We simulate the following background processes (with the generated number of events in 
parentheses): ZZ (10^), WZ (10*^), WW (10^), tt (10^), Z+jets (8-10^) and VT+jets (8- 10*^). 

3.1 Backgrounds to the trilepton channel 

We start with the W Z background, which is known to be the major source of background 
for the 3L channel. The total WZ cross section at Run II will be ~ 2.6 pb. Folding in the 
branching ratios of W and Z to leptons, we get a 3L WZ background cross section production of 
46 fb. It has a reducible and an irreducible component. The irreducible component (~ 3 fb) is 
due to Z ^ T^T^ t^i' decays. The invariant mass of the resulting lepton pair from the tau 
decays is usually far from the Z-mass, in a region which is typical of the signal. Hence there is 
no obvious cut which can substantially reduce this part of the W Z background without at the 
same time reducing signal-to-noise. On the other hand, the remaining background (~ 43 fb) 
is reducible, since it arises from Z t^t~ decays. In this case the invariant mass m^+^- of the 
resulting lepton pair is equal to p| = {pz)^{Pz)fi, where {pz)^i is the 4-momentum of the parent 
boson. Most of the time the Z is produced nearly on-shell, p| ~ M|. Hence, the invariant mass 
cut \m£+£- — Mz\ > 10 GeV is very efficient in removing this source of background. However, 
the parent can also be off-shell, due to either the Z-width, or to WZ/W'j interference. In 
PYTHIA, where only the first effect is modeled, the lepton pair invariant mass distribution 
follows a Breit-Wigner shape. We find that roughly 10% of the 3L background events pass 
the dilepton invariant mass cut, thus bringing the reducible background cross section down 
to about 4.3 fb. This is almost a factor of two larger than the corresponding irreducible 
background cross section (compare to 2.6 fb). Since ISAJET does not incorporate either Z- 
tail effect, we find essentially no reducible background cross section from ISAJET after the 
dilepton invariant mass cut is applied^. This difference between ISAJET and PYTHIA largely 



accounts for the discrepancy in the backgrounds found in Refs. [10, 12 1 and Refs. M O, M. 



■^Energy smearing in the detector simulation produces a very small background which survives the 10 GeV 
Z-mass window cut. 
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Figure 1: The invariant mass distribution of any pair opposite sign, same flavor 
leptons for the signal events (with Mq = 700 GeV, M1/2 = 160 GeV, tan/3 = 5) 



and the PYTHIA W Z background. We impose a set of cuts from Ref. [14|: pt{^) > 
{11,7,5} GeV, central lepton with px > 11 GeV and |r/| < 1.0, > 25 GeV and 
\m£+£- — Mz\ > 10 GeV. Each histogram is normahzed to its cross section. 



In what follows we use PYTHIA for our background estimate. 

We illustrate the above discussion in Fig. |T] where we show the invariant mass distribution of 
any pair of opposite sign, same flavor leptons for both signal and WZ background. The signal 
point has Mq = 700 GeV, M1/2 = 160 GeV and tan/? = 5, which results in m-± ~ m^o ~ 122 
GeV. The leptons are required to pass the set of cuts from Ref. [0]: Pt(^) > {H, 7, 5} GeV, 
central lepton with px > 11 GeV and |?7| < 1.0, J^t > 25 GeV and Irrii+i- — Mz\ > 10 GeV. 
The histograms are normalized to the respective cross section. 

Even with the modeling of the Z-width effect, we caution that the WZ simulation in 
PYTHIA is still not realistic, since the W'j* contribution is neglected. It results in a peak at 
low invariant mass, and the resulting distribution is markedly different from the result shown 
in Fig. |l|. We anticipate that it will be necessary to cut away all events on the low-end of the 
dilepton invariant mass distribution. We therefore always apply the cut \Tni+g- \ > 11 GeV for 
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all three channels. This cut also helps in eliminating background lepton pairs from Drell-Yan, 
as well as J/ip and T decays. 



Since recent trilepton analyses of the Tevatron reach |T3|, |I^ use ISAJET for the simula- 
tion, they underestimate the trilepton background. As a result, we find a significantly reduced 
Tevatron reach. However, in some cases we can employ an invariant mass cut which will sub- 
stantially reduce the WZ background with no significant loss of signal. The signal distribution 
in Fig. 1^ has a sharp kinematic cut-off at around fn^o — m^o ~ m-±/2 ~ 60 GeV, and we 
can exploit this feature to increase signal-to-noise. Indeed, by applying the more stringent cut 
mi+i- > m-±/2 ~ 60 GeV, we can eliminate most of the off-shell Z events, at almost no cost 
to signal. This is why in addition to standard Z-mass window cuts we consider four dilepton 
mass cuts which eliminate all events above a given invariant mass value. Also notice that 
only a very small fraction of signal events have invariant dilepton masses between and 11 
GeV, so that the low-end invariant mass cut rrii+i- > 11 GeV is quite efficient in improving 
S/ Our discussion of the WZ background can be similarly applied to the less serious, but 
nevertheless non-negligible ZZ background. 

The second largest background to the 3L channel is from dilepton ti events, where there 
happens to be a third isolated lepton from a 6-jet. This background is most easily suppressed 
by a jet veto. 

Finally, the remaining 3L background one should worry about is Z-iet production. In 
this case part of the background is due to events where a jet fakes a lepton. Monte Carlo 
simulations, especially with a simplified detector simulation like SHW, cannot give a reliable 
estimate of this background. In order to estimate the fake rate one has to understand the 
details of the detector response as well as the jet fragmentation. Only with Run II data will 
one be able to obtain a good estimate. For our study we follow a procedure which makes use 



of Run I data. It was used in Ref. to study the W + jets background to the 2L channel 
(see the next subsection). 



9 



3.2 Backgrounds to the like-sign dilepton channel 

We now discuss the backgrounds to the 2L channel. Ref. |2C] observed that it can be advanta- 



geous to not require the odd-sign lepton in the 3L events due to the gain in signal acceptance. 

At the same time, events with two like-sign leptons are still quite rare at the Tevatron, so the 

2L channel was suggested as a possible alternative to 3L for SUSY searches in Run I. However, 

the rates for the relevant diboson backgrounds were somewhat underestimated, since ISAJET 

was used for the simulation. Here we are interested in determining whether the 2L channel 

will be useful in the larger background environment of Run II. 

We first do a back-of-the-envelope comparison of the WZ backgrounds for the 3L and 2L 

channels. After not requiring the odd-sign lepton, one is left with the choice of vetoing that 

lepton. In the case of a veto, we find for the relative size of the two backgrounds 

aiyz(2L) 2£,(l-£0 +0.35£, [0.65 + 0.35(1 I - Si 

— ~ 1 1 1 

(Twzi'il^) 2efEz + 0.35%2 siEz ' ^ ^ 

where ei is the acceptance for leptons coming directly from W or Z decays, Sr is the acceptance 

for the (usually softer) leptons coming from leptonic tau decays, and ez is the efficiency of the 

Z-window invariant mass cut: \mg+g- — Mz\ > 10 GeV (we have neglected its effect on the 

irreducible background component). We find from Monte Carlo that the typical Run II values 

for these efficiencies in WZ production are ei ~ 0.60, Sr ~ 0.46 and ez ~ 0.10. Plugging into 

Eq. (|l|), we find for the ratio 6.3, which agrees reasonably well with the result 5.8 from our full 

Monte Carlo simulation. For the signal point shown in Fig. |l] we find ei = 0.62 and ez = 0.99 

and the corresponding ratio is 

asignal{'2L) ^ ^ ~ ^ 
(^signal 

(3L) eiez 

This reveals that vetoing the third lepton is definitely not a good idea. In comparison to the 
3L channel, the signal goes down, while the major background component is increased almost 
6 times! 

We therefore only consider the inclusive 2L channel, where we do not have any requirements 



on the third (odd-sign) lepton, just as in Ref. |2^. In this case, the signal acceptance is 
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definitely increased. Unfortunately, the corresponding increase in the background is even 
larger than before: 

CTyy^ (2L) ^ 2£;(1 - £; + £;g^) + 0.35£^ 

aiyz(3L) " 2^25^ + 0.35%2 ~ " ^ ^ 

for the typical values of the efficiencies. We can see immediately that the 2L channel can 
compete with the 3L on the basis of S/ \/B only if the lepton acceptance for the signal is less 
than 1/V7.3 ~ 37%. However, for typical values of the SUSY model parameters the lepton 
acceptance is much higher. 

To make matters worse, the 2L channel suffers from a potentially large new source of 
background: ly+jet production where the jet fakes a lepton. Although the rate for a jet 
faking a lepton is quite small, on the order of 10~^, the large ly+jet cross section results in 
a major background for the 2L channel. As we mentioned earlier, the best way to estimate 
this background is from data, since Monte Carlo simulations are not reliable for fakes. In our 
analysis we shall follow the procedure of Ref. [Q, where the rate for observing an isolated 
track which would otherwise pass the lepton cuts was measured in the Run I Z+jet event 
sample. This rate was then multiplied by the probability that, given an isolated track, it 
would fake a lepton. This probability was measured in Run I minimum bias events to be 
~ 1.5%, independent of px P3. In our study we first simulate with Monte Carlo the pr 
distribution of isolated tracks in W and Z production, which is shown in Fig. ^ In this figure 
we plot the isolated track pt distribution corresponding to the 2L background, i.e. we combine 
a real lepton with px > H GeV and |?7| < 2 with a same sign isolated track with \ri\ < 2. 
Hence, the 2L background cross section is obtained by multiplying the cross section from Fig. |^ 
by the probability that an isolated track will fake a lepton. We see that the background from 
fakes falls extremely fast with px, so a larger px requirement will substantially suppress it. 

We normalize the isolated track rate to data. Using the measured 1.5% fake rate per 
isolated track, we find 1.5 fb of cross section when running the simulation at ^/s = 1800 GeV 



and using the set of cuts from Ref. |20[. This is half the cross section found in Ref. pO|. Hence, 
to match the data to PYTHIA/SHW we need to double the isolated track rate obtained from 
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Figure 2: pT distribution of isolated tracks in (a) M/^+jet production and (b) Z+jet 
production. The isolated tracks have I < 2 GeV and are outside the AR = 0.7 
cone around any jet. The events in the distributions potentially contribute to the 
2L background: they have one real lepton with px > 11 and |?7| < 2 and one same 
sign isolated track with |r/| < 2. 



Monte Carlo. 

3.3 Backgrounds to the dilepton plus tau jet channel 



The largest background to the 2L1T channel is Drell-Yan [12], where the tau jet is a fake. As 



it turns out, SHW does quite a good job in simulating the fake tau rate so we can safely 
rely on the Monte Carlo for this background. We shall not further discuss the backgrounds to 



the 2L1T channel; we refer the interested reader to Refs. [Il2l p8[. We find a marginal increase 



in the background rate here, due to the larger i] coverage used in this analysis. 
3.4 Summary and discussion 

In Table |] we summarize our results for the different backgrounds to the three channels. 
The results are presented for a standard choice of cuts in each case. For the 3L channel. 
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Process 


3L 


|14| 


2L 


20 


2L1T 








E = 2 TeV 

m\<'2 


E = 1.8 TeV 

m\<i 




ZZ 


0.21 ± 0.01 


0.04 


1.836 ± 0.006 


0.1 ± 0.01 


0.372 ± 0.003 


wz 


1.39 ± 0.01 


0.40 


8.79 ± 0.03 


1.1 ± 0.02 


1.36 ± 0.01 


WW 


0.009 ± 0.003 


— 


0.002 ± 0.001 


+ 0.02 


0.54 ± 0.02 


ti 


0.33 ± 0.01 


0.14 


0.21 ± 0.01 


+ 0.02 


1.64 ± 0.03 


Z+jet 


0.13 ± 0.01 




3.58 ± 0.05 


1.1 ± 0.1 


11.3 ± 0.6 


W+jet 






10.5 ± 0.2 


3.0 




total 


2.07 ± 0.02 


0.58 


24.9 ± 0.2 


5.6 


15.2 ± 0.6 



Table 2: Background cross sections after cuts (in fb) for the three channels, 
each with a set of cuts described in the text. All errors are statistical. We 



also list the 3L background found in Ref. |I^ and the 2L background found 
in Ref. m. 



we pick the set of cuts from Ref. [|I^: Pt(^) > {11,7,5} GeV, central lepton with px > 11 
GeV and |r/| < 1.0, > 25 GeV, \mg+g- - Mz\> 10 GeV, and no |Ay?| or JV cuts. For the 



2L channel we use the following cuts from Ref. ||2^: Pt(^) > {H; 11} GeV, — Mz\ > 10 

GeV, and no l^y^l or JV cuts. Here we require both leptons to have \r]\ < 2.0, whereas 



Ref. requires \ri\ < 1. Finally, for the 2L1T channel we use the cuts Pt{t) > 15 GeV, 
\ri{T)\ < 1.5, prii) > {8,5} GeV, \ri{i)\ < 2.0, \me+i^ - Mz\ > 10 GeV, > 20 GeV, and no 
JV, from Ref. |T^. All errors in the table are statistical. 



Comparing to Refs. [p!4| , pO] we see a significant increase in the WZ and ZZ backgrounds 
to both the 3L and 2L channels. In the 2L case this is in part due to the larger rj coverage that 
we use, as well as the higher center-of-mass energy. The remaining difference is due to the lack 
of Z-tail effects in the ISAJET simulation. The WZ background is almost as large as ly+jet, 
which is not surprising based on the estimate in Eq. (compare the WZ backgrounds to 3L 
versus 2L). We remind the reader that our 2L W^+jet background has been normalized to the 
2L py+jet background from Ref. |20|, so the difference in the 14^+jet background seen in the 
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table is due solely to the different tracking coverage and center of mass energy. 

In the 3L case, we find significantly larger backgrounds than Ref. In fact, we surmise 
that the lepton efficiency is smaller in SHW than in the ISAJET detector simulation, and if 
one were to take this into account the background differences would be even larger. Part of 
the difference in the WZ and ZZ backgrounds is due to the fact that PYTHIA gives a ~ 15% 
larger diboson cross section than ISAJET. However, the largest part of the difference in the 
diboson background rates is due to the Z-width. Notice that with the fake rate procedure 
discussed in Sec. |3.2| we are able to obtain an estimate of the Z+jet trilepton background 
(where the jet fakes a lepton). This background has not been taken into account in previous 
studies. 

We do not trust the Monte Carlo simulation to provide a reliable estimate for the TV+jet 
and Z+jet backgrounds where the jet gives rise to a real isolated lepton. We expect these 
backgrounds to be small, and we ignore them. 

In the next section, we present results for the Tevatron reach in the three channels in the 
minimal gravity-mediated SUSY breaking model. In our simulations we use PYTHIA for the 
background estimates and ISAJET for the signal. 

4 Discovery reach for gravity- mediated models 

We next discuss the discovery potential of the upgraded Tevatron in the mSUGRA model. 
There are various SUSY production processes which give rise to the lepton signals under consid- 
eration. The dominant source of signal in most regions of parameter space is X1X2 production. 
Other chargino/neutralino and slepton production processes also contribute. Typically these 
processes constitute a small fraction of the total signal cross section, but in some regions of 
parameter space (e.g. large tan (3, small Mq) they can dominate. We ignore the possibility of 
small contributions from squark and gluino production processes. 

We first review which parts of the mSUGRA parameter space are accessible in Run II via 
the 3L signature. The X1X2 production cross section depends primarily on the chargino mass. 
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and scales roughly as M^.^ . Therefore, at the Tevatron we can only explore regions with 
small Ml/2, where the cross section is large. The Xi')<^ production proceeds predominantly 
via s-channel VT-boson exchange. There is also destructive interference from t-channel squark 
exchange. As a result, the Xi^^ cross section is slightly enhanced in case of heavy squarks 
(i.e. at large values of Mq). 

At large values of Mq (Mq > 700 GeV) the sleptons are heavy and the gauginos decay 
to three-body final states dominantly via real (at large M1/2) or virtual (small M1/2) W- or 
Z-boson exchange. In this case the reach is determined solely by the signal production cross 
section. At smaller values of Mq the off-shell slepton mediated decays destructively interfere 
with the gauge mediated decays. At small values of Mq (Mq ^ M1/2) sleptons become lighter 
than the Xi I X2 ^"^^ two-body decays to lepton final states open up. This enhances the leptonic 
branching ratios of the gauginos and improves the Tevatron reach. We illustrate these features 
in Fig. ^ where we show the branching ratios|of a chargino-neutralino pair into the 3L, 2L and 
2L1T channels at an mSUGRA model point with tan/5 = 5, /i > 0, = and (a) M1/2 = 175 
GeV or (b) M1/2 = 250 GeV. We plot versus m^^ (bottom axis) or Mq (top axis). Notice that 
because of the rather small value of tan /5, all slepton fiavors are practically degenerate. 

In Fig. ^(a) the chargino mass m-± is only about 120 GeV, so the two-body decay to 
the W^-boson is closed. If Mq > 92 GeV the chargino decays are three-body and the 2L and 
3L channels have a larger branching ratio than the 2L1T, roughly by a factor of two. The 
dip in the leptonic branching ratios near m^^ ~ 260 GeV is due to destructive interference 
between the Z and ^-mediated graphs In the region Mq < 92 GeV the two-body decays 
to sleptons are open and quickly become dominant. But notice that while X2 can decay to all 
slepton fiavors through its bino component, xf decays dominantly to a stau, since the chargino 
couplings to the right-handed sleptons are proportional to the corresponding lepton masses. 
Hence, the branching fraction of the 2L1T channel dominates the region 40 ^ Mq ^ 90 GeV. 



§When we use the term 'branching ratios' in relation to a pair of particles, we mean a smii of products 
of branching ratios. For example, if the two-body decays of X\ ^-^d X2 ^-''^ closed, BR(x^X2 ~' 3i) = 
[BR(x^ ^ Xl^^yd + BR(x^ ~- x\T^Vr ^ x?^+^^i^rz^r)] [BR(x^ ^ x\t^n + BR(x§ ^ x\t^t- ^x? 

t^l Vi ViVr Vt)]- 
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Figure 3: Branching ratios of a chargino-neutralino pair into the 3L (sohd), 2L 
(dash) and 2L1T (dot dot dash) channels versus the lightest stau mass (bottom axis), 
or alternatively, versus Mq (top axis), with mSUGRA model parameters tan/3 = 5, 
/i > 0, Ao = and (a) M1/2 = 175 GeV or (b) M1/2 = 250 GeV. The arrows indicate 
the chargino threshold m-± = mf, . For the range of Mq values shown the chargino 
mass varies from 118 to 139 GeV in (a), and from 187 to 201 GeV in (b). 



For Mq < 40 GeV the 2L1T branching ratio rapidly drops while the 2L and 3L branching 
ratios increase and become dominant again, due to two-body decays of Xi to sneutrinos. 

In Fig. 11(b) the mass of Xi is near 190 GeV and it is always above the VT-boson threshold. 
If Mo is greater than about 400 GeV the X2 is above the Z-boson threshold, and in this region 
the branching ratios for the three signatures are determined by the W- and Z-boson branching 
ratios to leptons. In the region Mq < 150 GeV the sleptons are lighter than the chargino, 
and the two-body decay modes to sleptons greatly enhance the leptonic branching fractions. 
The 2L1T channel does not become as dominant because the chargino decays to quarks via 
an on-shell VT-boson. The decrease in the leptonic branching ratios below Mq ~ 70 GeV is as 
before due to two-body decays to sneutrinos. 

In Fig. ^ we show similar plots of the branching fractions, but for tan /5 = 35. Many of the 
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Figure 4: The same as Fig. ^, but for tan /3 = 35. This time the chargino mass 
varies from 126 to 140 GeV (in (a)), and from 192 to 201 GeV (in (b)). 



features seen in Fig. ^ are present here as well. For example, in Fig. ^(a) we see the broad 
region (150 ^ m^^ < 400 GeV) where the destructive Z-i interference severely diminishes the 
branching ratios. And again we observe the X2 ~^ Xi^ threshold in Fig. ^(b), near nif^ = 260 
GeV. However, in this case the tau slepton is significantly lighter than the first two generation 
sleptons. Hence, the branching ratio to the three tau final state quickly approaches 100% below 
the stau threshold, and the 2L1T branching ratio is large. We also see that the 2L channel is 
competitive at small values of Mq, and is preferred over 3L on the basis of branching ratio by 
a factor of 2.8. 

To summarize, we see that depending on the values of the mSUGRA parameters, any of 
the three channels offers some promise to be observed in Run II. In the rest of this section, 
we shall do a comparative study of the three channels, accounting for all relevant background 
processes and using a realistic detector simulation. We shall not only update the existing 
3L analyses with improved estimates of the WZ, ZZ and Drell-Yan background rates, but 
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foremost we are interested in evaluating the different prospects each channel can offer. For 
example, we would like to see whether the 2L channel offers reach beyond the 3L channel, or 
can cover regions where the 3L channel is suppressed. Also, we want to determine the region 
of parameter space where the 2L1T channel can offer an independent check of a signal in one 
of the other channels. 

The cut optimization procedure is an important ingredient in our analysis. We show results 
for the reach with the optimal set of cuts, determined independently at each point in SUSY 
parameter space. The optimal set of cuts maximizes S/ \/B. We require the observation of at 
least 5 signal events, and present our results as Sa exclusion contours in the Mq — ^1/2 plane, 
for two representative values of tan/? - 5 and 35. We fix /i > and = 0. 

In Fig. ^ we show the Tevatron reach in the 3L channel with a standard set of soft cuts 
(dashed lines, for large Mq only to prevent crowding), as well as with the optimal set of 
cuts (solid). We show the expected reach for 2, 10 and 30 fb~^ total integrated luminosity. 
The cross-hatched region is excluded by current limits on the superpartner masses and the 
dot-dashed lines indicate the projected LEP-II reach for the chargino and the lightest Higgs 
massPj. In Fig. |^(a) the left (right) dotted line marks the W/xf i^i/xf) mass threshold, while 
in Fig. 1(b) the left (right) dotted line marks the eu/xf (^i/x^) mass threshold. 

Comparing to the results of Refs. [|l^, |l^, we see the Tevatron reach indicated in Fig. |^ 



is greatly reduced, due to the larger backgrounds found in our analysis. Even with optimized 
cuts our results show a much smaller observable region. For example, using the set of cuts of 



Ref. [|14| at tan j3 = 35 and Mq = 1 TeV, our results indicate that the region bounded by the 
2 fb^^ 3-cr contour extends to M1/2 = 123 GeV (or 136 GeV with optimal cuts), whereas in 
Ref. the corresponding 5-0" region extends to 180 GeV. Similarly, their 30 fb~^ 3-0" contour 



extends to M1/2 = 250 GeV, while ours extends to M1/2 = 186 GeV (198 GeV with optimal 
cuts). 

^It should be kept in mind that at small tan/? the Higgs mass is a sensitive fmiction of tan/3. For example, 
if the Higgs-boson is not discovered at LEP-II, tan/3 = 3 will be excluded. The reach contours are much less 
sensitive to tan f3. 



18 




200 400 600 800 200 400 600 800 1000 

Mo (GeV) Mo (GeV) 

Figure 5: Tevatron reach in the 3L channel for mSUGRA models with > 0, 
^0 = 0, and (a) tan/? = 5 or (b) tan/5 = 35. We show the reach with both 
a standard set of soft cuts [|l4| (dashed, for large Mq only), as well as with the 
optimal set of cuts (solid) (see text). The reach is shown for 30 fb~^, 10 fb^^ and 
2 fb~^ total integrated luminosity (from top to bottom). The cross-hatched region 
is excluded by current limits on the superpartner masses. The dot-dashed lines 
correspond to the projected LEP-II reach for the chargino and the lightest Higgs 
masses. In Fig. (a) the left dotted line shows where niu^ = m-± and the right dotted 
line indicates mf^ = m~± . In Fig. (b) the dotted lines show where me„ = m~± (left) 

Xi Xi 

and mf. = m~± (right). 

Xi 

As expected, there is respectable reach beyond LEP-II at small values of Mq and tan/3, 
where the chargino and neutralino 2-body decays to sleptons are open. As expected from 
Fig. 0(a), the Tevatron has no sensitivity beyond LEP-II in the region 200 < Mq ^ 400 GeV. 
For large values of Mq, where the leptonic decays of the gauginos are three-body, we find some 
sensitivity for both values of tan (3. In this region there is a clear benefit to using optimized 
cuts. At Run II (2 fb~^), with the default set of cuts, only the region with small tan/3, small 
Mq, and small M1/2 can be explored beyond LEP-II. With optimal cuts, we see at large tan j3 a 
non-negligible region can be excluded beyond LEP-II. Looking beyond Run II, we notice that 
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Figure 6: The difference AL (in fb~^) between tlie required total integrated lumi- 
nosity for the optimal set of cuts and the default set of cuts, for the 3L signal. The 
three (solid) contours correspond to (from top to bottom) AL = 30, 10 and 2 fb~^. 
The dotted lines indicate the optimal-cut 30 fb~^ reach contours from Fig. ||. 



at TeV33 (30 fb^^) optimization can prove equivalent to doubling and sometimes even tripling 
the total integrated luminosity! In Fig. |^ we show the difference in the required luminosity 
when using the fixed set of cuts from |]T^ and the optimized cuts. As a guideline, we also 
show the 30 fb~^ optimal-cut contours (dotted lines) from Fig. ^. 

It is instructive to examine the optimized sets of cuts. In Fig. |^ (Fig. ^ we show the 
optimal sets of cuts for the 3L channel in the Mq, M1/2 plane, for tan/5 = 5 (tan/? = 35). We 
use the following notation to describe the set of cuts at each point. The central symbol is the 
number of the pt cut according to Table 0. The left superscript indicates which central lepton 
7] cut was chosen, either |?7| < 1.0 (labeled "10") or \ri\ < 1.5 ("15"). A left subscript denotes 
the type of \Aip\ cut: for the two highest pt opposite sign, same flavor leptons the \Aip\ < 2.5 
cut is indicated by "2", the \ A(p\ < 3.0 cut by "3", and no symbol indicates no \ Aip\ cut. The 
cut \A(p\ < 2.5 for any pair of opposite sign, same flavor leptons is indicated by "a". The 
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Figure 7: The optimal sets of 3L cuts in the Mq, M1/2 plane, for tan/3 = 5. The 
key indicates which symbols correspond to which cuts (central lepton rapidity, Ipx^ 
invariant mass, A99, and Pt cuts) (see text). The dotted line indicates the optimal- 
cut 30 fb~^ reach contours from Fig. ^. 

right superscript shows the $t cut: $t > {15,20,25} GeV ("15", "20", "25"), or no cut (no 
symbol). A right subscript denotes the dilepton invariant mass cut: \mi+i- — Mz\ > {10, 15} 
GeV ("10", "15") or |m£+£-| < {50,60,70,80} GeV ("50", "60", "70", "80"). And finally, a tilde 
over the central symbol indicates that the luminosity limit came from requiring 5 signal events 
rather than 3a exclusion. 

The jet veto cut is not indicated on the figures. However, except for one point at large 
tan/3 the jet veto was never selected as an optimal cut. Indeed, the major 3L background 
events (from WZ) are just as likely to contain extra jets as the signal. 

There are several lessons to be drawn from Figs. |^ and ^. As expected, in those cases where 
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Figure 8: Same as Fig. 0, with tan /? = 35. 



the signal is strong and dominates over the background (typically for small M1/2 or small Mq), 
softer cuts are beneficial. Indeed, the majority of the points which can be discovered with 5 
signal events, have selected px cut "1" , which is the most lenient set of px cuts. For larger 
values of M1/2, where the background is more important, harder cuts on the lepton p^^s are 
preferred. In fact, we see that in the region of interest for TeV 33 at large Mq, the hardest pt 
cuts we consider (from Ref. [0, indicated by "5"), often work best. 

Figs. 1^ and |^ clearly indicate the advantage of the invariant mass cuts introduced in 
Section |^. Of all the points on the plots where the background is an issue, there are only a few 
at which a conventional Z-mass window cut is optimal. In all other cases it is advantageous 
to cut all events with invariant dilepton masses above a certain threshold. Notice also how 
the value of the threshold tends to increase with increasing M1/2. This is expected because 
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the sharp edge in the signal distribution is roughly at 0.4 ■ M1/2 (see Fig. p. 

We also see from the figures that a missing cut is preferred essentially everywhere in 
parameter space. However, a soft cut > 15 GeV) often gives the better reach. The 
latest trilepton analyses [|1^, have chosen to increase the cut from its nominal Run I 
value of 20 GeV to 25 GeV. Our results suggest that in the off-line analysis one is better off 
with a lower cut. However, more work is needed to conclusively determine what the best 
fix cut will be in the actual analysis. For example, triggering and energy mismeasurement 
issues will have to be carefully taken into account. 

Lastly, in the great majority of parameter space it is better not to require a ip cut. 

Looking back at Fig. |^(a), we observe that near the slepton-chargino mass thresholds the 
reach becomes diluted. This effect can easily be overlooked, since it only shows up very close 
to threshold. We find that it is entirely due to the suppressed signal acceptance. Indeed, 
although the branching ratio is increased immediately below threshold, the lepton resulting 
from the X2 ~^ ^^^^ decay tends to be very soft and it can fail the analysis cuts. 

We next discuss the prospects for the 2L channel. In Fig. |^ we show the 2L channel reach 
for the Tevatron. Due to the much larger background the reach in the 2L channel is not 
as great as in the 3L channel. There is one important exception, however - the 2L channel 
does not lose sensitivity near the chargino-charged slepton threshold. Indeed, because of the 
Majorana nature of the neutralinos, in 50% of the events the lost soft lepton in the X2 ~^ i^i^ 
decay is of the opposite sign as the charged lepton from the chargino decay. The remaining 
two hard leptons are then of the same sign and they are readily reconstructed. Therefore, 
near the i/xf mass threshold, the 2L channel may prove to be a valuable alternative to 3L. 

Fig. ^(b) reveals that at large tan (3 one starts to lose 3L sensitivity in the small Mq region, 
since the decays to tau final states dominate. In fact, for Mq ^ 300 GeV and tan/5 = 35, 
we find no 3L reach in Run II beyond LEP-II. Only with multiple years of running and 
collecting soft lepton events will the Tevatron be able to start improving on the LEP mSUGRA 
bounds. One can improve this situation by considering alternative signatures with tau jets 
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Figure 9: The same as Fig. ^, for tlie 2L channel. The dashed lines correspond to 
a set of cuts from Ref. [2C]. The 2, 10 and 30 fb~^ reach contours are plotted. 



| T^ . Of those, the 2L1T channel is singled out on the basis of both sensitivity and statistical 
importance. 



In Fig. [ly we show the Tevatron reach in the 2L1T channel. In this case we compare the 
optimal reach to a set of cuts from Ref. |]12| (listed in Sec. [3^ ). We find very little difference 
between this fixed set of cuts and the optimal set, so we do not show the fixed cut lines in the 
figure. 

The 2L1T channel has no reach in the large Mq region. However, when the two-body 
decays to staus open up the 2L1T branching fraction is larger than the other two channels, 
leading to some sensitivity. While the region accessible in the 2L1T channel at small tan/3 is 
not competitive with the 3L and 2L reach, it can improve the statistical significance in case of 
exclusion, or it can serve as an important confirmation and provide unique information about 
the model parameters in case of discovery. This channel offers the greatest reach at large tan (3 
(together with the related much cleaner signature of two like-sign leptons plus a tau jet [0) 
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Figure 10: The same as Fig. |5|, for the 2L1T channel. In this case the optimal cuts 
yield very little improvement relative to a default set of cuts from Ref. ||l^ (listed 
in Sec. |3.4|) , so the contours corresponding to the default cuts are not shown. 

in the small Mq region. 

5 Conclusions 

In this paper we studied three of the cleanest and most promising channels for SUSY discovery 
at the Tevatron in Run II. We revisited the trilepton and like-sign dilepton analyses, improving 
them in several key aspects. For example, we used a more realistic simulation of the major 
backgrounds. We found larger backgrounds than previous analyses. We used a procedure 
relying on Run I data to estimate backgrounds involving fake leptons. And we introduced 
an invariant mass cut which took advantage of a sharp edge in the signal dilepton invariant 
mass distribution. This cut was generally more effective in increasing signal-to-noise than the 
standard invariant mass cuts. Also, we varied the cuts at each point in the supersymmetric 
parameter space, and determined at each point the set of cuts which yields the largest reach. 
We found that this cut optimization can significantly enhance the Tevatron reach. Lastly, we 
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analyzed the reach of the 2L1T channel. 

If nature is supersymmetric at low energies, and the superpartners (in particular the gaug- 
inos) are light, there is a good chance that the Tevatron will discover them in its upcoming 
runs. However, the 3L signature has limited reach, and we can improve the reach by optimizing 
cuts and considering alternative clean signatures. 
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